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ABSTRACT 

We report the results of a 944-epoch survey for transient sources with archival 
data from the Very Large Array spanning 22 years with a typical epoch sep- 
aration of 7 days. Observations were obtained at 5 or 8.4 GHz for a single 
field of view with a full-width at half-maximum of 8.6' and 5.1', respectively, 
and achieved a typical point-source detection threshold at the beam center of 
^ 300 /iJy per epoch. The angular resolution ranged from ~ 1" to 15". Ten 
transient sources were detected with a significance threshold such that only one 
false positive would be expected. Of these transients, eight were detected in 
only a single epoch. Two transients were too faint to be detected in individual 
epochs but were detected in two-month averages. None of the ten transients 
was detected in longer-term averages or associated with persistent emission in 
the deep image produced from the combination of all epochs. The cumulative 
rate for the short timescale radio transients above 370 /iJy at 5 and 8.4 GHz is 
0.07 ;^ -R ;^ 40 deg^^yr^^, where the uncertainty is due to the unknown duration 
of the transients, 20 min ^ tchar ;^ 7 d. A two-epoch survey for transients will 
detect 1.5 ±0.4 transient per square degrees above a flux density of 370 /iJy. One 
transient is located ~ 3 kpc in projection from the nucleus of a spiral galaxy 
at a redshift z = 0.040. Based on the duration of the transient, its luminosity 
(L ^ 2 X 10'^^ erg s~^), and its location in the galaxy, we suggest that it may 
be similar to the peculiar Type Ib/c radio supernova SN 1998bw associated with 
GRB 980428. The implied Type Ib/c rate is consistent to an order of magnitude 
with rates determined optically. A second transient is associated with a blue 
galaxy at z = 0.249 and may also be a luminous radio supernova or gamma-ray 
burst afterglow. Two other transient sources are possibly associated with the 
outer parts of faint, optically detected galaxies. The remaining six transients 
have no counterparts in the optical or infrared to limiting magnitudes of i? ~ 27 
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and Kg ~ 18 and show no faint persistent radio flux. For the eight transients 
without clear associations, known source classes of radio transients including ra- 
dio supernovae, gamma-ray burst afterglows, active galactic nuclei, and stellar 
flares do not provide natural fits to the observed parameters. The hosts and 
progenitors of these transients are unknown. 

Subject headings: radio continuum: general — radio continuum: stars — super- 
novae: general — gamma rays: bursts — surveys 



1. Introduction 



Radio transient sources have primarily been studied th rough follow-up observations af 



ter discovery at optical, X-ray, or 7-ray wavelengt 
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Bower et al. 2003). This 



is principally due to the high cost of observing time t o survey large areas of sky to sig- 
nificant depth at radio wavelengths (ICordes et al.ll2004l ). Moreover, existing blind surveys 
have typically been performed at low frequencies, where survey time is shorter but the ef- 
fect of synchrotron self-absorption may hide certain source classes. The available parameter 
space for radio transient surveys is extensive: transients have been detected at, and are 
predicted for all, radio wavelengths ; timescales of transient and variable behavior range from 
nanoseconds ( Hankins et al.l I2OO3I ) to the longest timescales probed ( Hughes et al. 1992); 
and transie nts may originat e from nearly all astro physical environme nts including the so 
l ar system (Katz et al. 200c 



( IZhao et al. 



1992 



star- fo rming regions ( Bower et al. 2003), th e Galactic Center 



Bower et al.ll2005l ). and other galaxies (IFrail et al.lll997l ). 



Nevertheless, extensive surveys have been carried out, exploring different volumes of 
the parameter space, and, in some cases, discovering new phenomena. Roughly, one can 
separate transients surveys into two classes: (a) burst searches that probe timescales of less 
than ~ 1 second, are often performed with low angular resolution instruments, and are most 
often performed at low frequencies; and (b) imaging surveys conducted with interferometers 
which typically probe timescales of tens of seconds and long er. Examples of burst searches 
include STARE, a 611-MHz all-sky survey (IKatz et al.ll2003l ) sensitive to timescales of 0.1 s 
to a few minutes; a 843-MHz survey with the Mol onglo Observator y synthesis Telescope 
(MOST) sensitive to timescales of 1 /is to 800 ms (lAmy et al.lll989l ): and the Parkes 1.4- 
GHz multi-beam pulsar survey. STARE used dipole antennas with a sensitivity threshold 
of 27 kJy simultaneously over a significant fraction of the sky and found no extra-solar 
transients on timescales of 0.125 s to a few minutes in 18 months of observations. The 
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MOST survey determined an upper limit to the transient rate of 5 events s~^ deg~^ for 
10 ms events at a flux density limit of 1 Jy. The Parkes multi-beam survey has been used 
for extensive searches for periodic emission as well as for pulsed emission. The latter search 
led to the discov ery of rotating radio tra nsients (RRATs), which appear to be pulsars with 
erratic emission (IMcLaughlin et al.ll2006l ). 



Imaging surveys for transients have focused on individual fields and on large-scale sur 
veys. Recently, a blind survey of the Galactic center at 330 MHz has discover ed an unusual 



radio t ransient that displays quasi-periodic emission on a timescale of 70 min (IHyman et al. 
20051 ). ICarilli et al.l (120031 ) used ~ 30 hr of observation of the Lockman hole to character- 
ize variability on a timescale of 19 d and 17 mon at 1.4 GHz. They found no transients 
above 100 /xJy and conclude t hat the transient density on these timescales is fewer than 18 
per square degree. Similarly, iFrail et al.l (120031 ) identified 4 highly variable sources but no 
transient sources at 5 and 8.4 GHz using VLA observations of fields in which gamma-ray 



bursts are present. They set an upper limit to the variable source po 



sub-mJy sensitivity, comparable to the limit set by lCarilli et al.l (120031 ). Recently, drift scan 



Dulation of 6 deg ^ at 



interferometric observations have identified several transients of unknown origin, including 
a 3-Jy transient with a duration of 72 hr and three > 1 Jy transie nts of i ^ 1 d duration 
jNiinuma et al.ll2007l : iKunivoshi et al.ll2007l : iMatsumura et al.l 120071 ). 



The most extensive imaging se arch for transients derives from a comp arison of the 
FIRST and NVSS 1.4 GHz catalogs Jlevinson et al.l I2OO2I : ICal- Yam et al.lbooeh . The search 
covered 2400 square degrees at a flux density threshold of 6 mJy. This search generated a 
number of transient candidates. Of these, one source was identified as a radio supernova in a 
nearby galaxy and another was seen only once in the radio and could not be associated 
with any known optical source. Transient identification in this survey suffers from the 
methodological problem of the extremely mismatched resolutions (FIRST: 5", NVSS: 45") 
of the two surveys. 

In addition to radio supernovae (RSNe), a number of highly variable and transient 
sources are expected to be detectable in significant numbers at radio wavelengths. One 
example is gamma-ray burst (GRB) afterglows associated with GRBs which are b ear ned in a 
direc tion away from Earth, and therefore not observable at high-energy wavelengths (IRhoads 
19971 ). These so-called "orphan gamma-ray burst afterglows" (OGRBAs) are expected to 
outnumber GRBs by a factor of 10 to 1000, the specific value depending on the characteristics 
of the relativistic outflow of the GRB and the surroundi ng medium. There ni ay be related 
phenon iena such as those associated with X-ray flashes (ISakamoto et al.l 120051) . supernovae 



(SNe) ( Berger et al. 2003), and short-duration GRBs ( Gehrels et al.l 20051 ) . Activity from 



stars and compact objects in the Galaxy may also be detected as transients (e.g. jBerger et al. 



-4- 



200ll : Bower et alJboOSl : lOsten et al.ll2006h . 



A range of active galactic nucleus (AGN) phenomena is also anticipated, including intra- 
day variability and other short-term variability driven by the e ffect of interstellar s cintillation 
on compact components of extragalactic radio sources (e.g., iLovell et al.l 120031 ). Intrinsic 
AGN variability is also expected. Sudden increases in the mass accretion rate following the 
tidal disruption of an or biting star, for instance, will pr oduce a transient soft X-ray flare and 
possibly radio emission (lReeslll988l : iGezari et al.ll2006l ). 



Significant numbers of RSNe, OGRBAs, and tidal flares are predicted for surveys that 
achieve a sensitivity ~ 0.1 mJy and cover an area of ~ 10 square degrees. These phenomena 
all produce self-absorbed synchrotron radiation and, therefore, are more likely to be observed 
at high radio frequencies. 

We describe here the results of a search for extragalactic radio transients using 5 and 8.4 
GHz observations from the Very Large Array data archives. The survey includes 944 epochs 
over 22 years consisting of observations of the same field. At a sensitivity threshold of 370 
/iJy, the survey has an effective survey area comparable to that of a survey consisting of 
two epochs of 10 square degrees each, making this thus far the largest-area sub-mJy imaging 
survey of radio transient phenomena. 



2. Observations and Analysis 

2.1. The Archival Survey 

Between 1983 and 2005, the VLA performed calibration and system check observations 
of the same blank field with a typical interval of ~ 7 d. Each epoch typically was 20 min 
in duration and consisted of 10 s integrations. The field is centered at a = 15^02^20. 53*, 
S = +78°16'14.905" (J2000). The field is out of the plane of the Galaxy (/ = 115°, b = 36°). 
626 epochs of observations were obtained at a frequency of 5 GHz between 1983 and 1999. 599 
epochs were obtained at a frequency of 8.4 GHz between 1989 and 2005. VLA observations 
of this field ceased in 2005. Between 1989 and 1999, there was overlap between some of the 
5 and 8.4 GHz observations. 281 of the 8.4 GHz epochs were simultaneously observed at 5 
GHz, giving a total of 944 independent epochs at either frequency. 

Data were all obtained in standard continuum mode with 50 MHz of bandwidth in each 
of two intermediate-frequency bands (with centers separated by 100 MHz) and two circular 
polarizations. We used a pipeline procedure for fiagging, calibrating, and imaging. Phase 
calibration was made based on brief observations of the compact source J1803+784. Since no 
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standard flux density calibrator was used in these observations, we set the amplitude scale 
of the observations by assuming the mean flux density for J1803+784 as measured in the 
University of Michigan Radio Astronomy Observatory dat abase over the same period (2.2 Jy 
and 2.8 Jy at 5 and 8.4 GHz, respectively AUer et aDllQod ). The assumption of constant flux 



density introduces less than 15% uncertainty in the flux-density scale for individual epochs. 
Images of the target field typically had a root-mean-square (rms) error in the flux density 
of 40 to 50 /iJy. Images from two successive epochs (19840613 and 19840620) are shown in 
Figured! Throughout the paper, we use the notation of YYYYMMDD to denote epoch date 
and refer to specific transients as RT YYYYMMDD, based on the epoch of detection. 

Images were obtained in all configurations of the Very Large Array. For the extended 
configurations, bandwidth smearing within the primary beam can become significant. We 
compensated for this effect by applying a Gaussian taper with a full-width at half-maximum 
(FWHM) of 150 kX to the {u, v) data before imaging, which reduces sensitivity and angular 
resolution but provides imaging to large radii without bandwidth distortions. Images were 
made using natural weighting to ac hieve maximal sensitivity and were deconvolved using 



the CLEAN algorithm (IClarkI Il980l ). We imaged a region slightly larger than the region 
containing a circle of radius equal to twice the half-power radius, which is 9.0' at 5 GHz and 
5.4' at 8.4 GHz. 

A deep image of the region at 5 and 8.4 GHz was also constructed through merging of 
all visibilities and imaging with a {u, v) taper of 150 kX using the MIRIAD software package. 
The resulting images had resolutions of ~ 5" and ~ 3", and an rms of 2.6 /xJy and 2.8 /zJy at 
5 and 8.4 GHz, respectively. The location of steady and transient radio sources on the deep 
image is shown in Figure [2l Fluxes in the deep image are typically lower than the average of 
detections since our high detection threshold and the low flux densities of the sources bias 
the mean. Further details of the deep images will be published separately. 



2.1.1. Transient Identification 

Sources were identified using AIPS task SAD (search and destroy), which identifies 
peaks in the image above a flux-density threshold and fits point sources. For the relatively 
uncluttered images that we have, SAD is efficient and accurate. We identified sources within 
twice the half-power radius and corrected their flux densities for primary beam losses. 

The detection threshold for each image was chosen such that the probability of a false 
detection (PFD) in the full imaged region was 10^^, or a single false source for 1000 identical 
epochs, assuming Gaussian statistics for noise in the image. Since the synthesized beam 
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changes with array configuration while the primary beam remains the same, the number of 
independent pixels in the field of view varied substantially. Detection thresholds were in the 
range ~ 5a to 6a, dependent on array configuration. The typical fiux density threshold at 
the center of the image was ~ 300 /iJy. For the number of epochs observed and the circular 
area of the region searched for sources, we have an expectation of ~ 1 false source in the 
entire survey. Inverting the sign of the intensity in each image and searching for sources 
with the same thresholds and methods revealed one significant source, consistent with our 
expectations of one or fewer false detections in the entire survey. 

Wider-field images (~ 40' at 5 GHz, ~ 27' at 8.4 GHz) were made for all epochs in 
which a transient was detected. A number of transients that were detected on the western 
and northern edges of the individual epoch image were determined to be aliased power from 
two bright radio sources (~ 50 mJy) at radii greater than two primary beam widths (NVSS 
J1457+7817, J1500+7827). One source above the detection threshold was rejected on the 
basis of a non-pointlike image and deep negative stripes in the image. 

Detected sources and their measured properties are listed in Tables [T] and [21 Table [T] 
identifies persistent sources above or near the detection threshold and not multiple-epoch 
transients. Six persistent sources at 5 GHz and 2 persistent sources at 8.4 GHz were detected 
in multiple epochs. The epochs of detection were distributed over the entire survey. Of these 
persistent sources, three radio sources with fiux densities > 500 /iJy were detected hundreds 
of times. One radio source (J150123-I- 781806) was detected with our pipeline software at 5 
GHz in 452 of 626 epochs, and at 8.4 GHz in 159 of 599 epochs (Figure E]); with a reduced 
source detection threshold of 3a, we find that this source was present in 602 epochs at 5 GHz 
and in 445 epochs at 8.4 GHz. The results for this source indicate the overall quality of the 
data. Variations in the fiux density of this source are caused by the changing sensitivity of 
the VLA in different configurations, variability in the amplitude calibrator J1803+784, and 
possibly intrinsic variability in the source. 

Seven sources at 5 GHz and 1 source at 8.4 GHz were detected in individual epochs 
(Table [2]). These sources were detected only once and are classified as transients. We plot 
the light curves of each transient for the year surrounding their detection (Figure H]). 

We searched for time variability for each transient by splitting the data into four-minute 
segments and comparing the fiux densities. The reduced fo^' the hypothesis of no vari- 
ability is less than 1.4 for all the transients, implying no strong evidence for variability. We 
also imaged each transient in Stokes V and found no evidence for circular polarization above 
~ 30%. Finally, we differenced images between the intermediate frequency bands and found 
no evidence for a large spectral index. 
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The radial distribution of sources is consistent with a real source population rather than 
imaging defects. Imaging defects would be uniformly distributed throughout the image, 
without regard for the decline in sensitivity due to the primary beam. Four of the eight 
transients are contained within the half-power radius. This number is inconsistent with the 
expectation of a uniform distribution within the same region, 2.0 ± 1.4. On the other hand, 
a source population with a power-law distribution of flux densities N{S) oc 5""^/^ in a fiux- 
density limited sample has an expectation of 3.6 ± 1.3 within the half-power radius, in better 
agreement with the observed number. 

Images were also constructed in a search for faint radiation from the transients and for 
new transients using all data from two-month, one-year, and one-decade time spans. We 
found two new transients in the two-month averages (Table [3]) and no transients in the year 
and decade images. Light curves for these sources are shown in Figure |H Given the smaller 
number of epochs, the expected number of false sources in the two-month average data is 
~ 0.1. Further searches on different time scales or with filters matched to light curve shapes 
are possible with the same data but are beyond the scope of this paper. 

Dual frequency observations were made only after 1989 and not at every epoch. Of 
the remaining 5 GHz transients, only 1 (RT 19920826) is located within the imaged region 
of the 8.4 GHz survey. Unfortunately, there were no 8.4 GHz observations for this epoch. 
Similarly, for RT 19970205 which was detected at 8.4 GHz, there was not a simultaneous 5 
GHz observation. 



2.2. Keck Imaging 



On 25 July 2006 (UT dates are used throughout this paper), we observed t he central 



porti on of the radio field with the Low- Resolution Imaging Spectrometer (LRIS; lOke et al. 



19951 ) on the Keck-I 10-m telescope. Eight of the radio transient positions (see Figure 6) 
were observed in two separate pointing centers (each covered by the LRIS field-of-view of 
6' X 7.8'). At each pointing center we obtained six dithered 10-min exposures. The instrument 
is outfitted with a beam splitter, allowing us to observe simultaneously in the g and R 
bands. The data were reduced following standard optical imaging reduction procedures. 
After fitting a world coordinate system (tied to t he USNO Bl . O) to each reduced frame, we 
made a mosaic for a given filter using SWARP (IBertin et al.l 120021 ). The final astrometric 
uncertainty relative to the International Celestial Reference System (ICRS) is 250 mas in 
each axis. 



On the night of the Keck/LRIS imaging, we observed the Landolt standard star field 
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PG1633+099 (iLandoltl Il992l ) in both filters at an airmass of sec z = 1.04. Tlie catalog 
magnitudes were converted to the g band from the \^-band magnitude and the B — V 
color using the prescription in Bilir et al. (2005) . From the observed photometry of four 
PG1633+099 stars, the zeropoint conversion from flux in both filters was determined in an 
0.6" aperture. The flux of 3(T-detected sources in both the R and g stacked image of the radio 
field were determined in the same-sized aperture. Accounting for the average extinction per 
unit airmass at Mauna Kea, we determine that the non-detections of some of the counterparts 
correspond to upper limits of 5'(lim) ^ 27.6 and i?(lim) ^ 26.55 mag. These limits assume 
an unresolved source and a color similar to those of the Landolt field stars. 



2.3. Keck Spectroscopy 



An optical spectrum of the counterpart of RT 19840613 was taken on 2005 Dec. 4 with 
the DEIMOS spectrograph (IFaber et al.l l2003l ) mounted on the Keck-II telescope, using a 
600 line mm~^ grating and the GG400 order-blocking filter. Spectra of the counterparts and 
nearby companions of RTs 19860115, 19870422, 19920826, and 19970528 were taken on 2006 
June 26 with LRIS using the 600/4000 grism, the 400/8500 grating, and the D560 dichroic. 



All spectral data were reduced using standard techniques (e.g.. lFoley et al.ll2003l ). Stan- 
dard CCD processing and spectrum extr action were co mpleted with IRAF0. The data were 
extracted with the optimal algorithm of iHornd (119861 ). We obtained the wavelength scale 
from low-order polynomial fits to calibration-lamp spectra. Small wavelength shifts were 
applied after cross-correlating night-sky lines extracted with the object to a template sky. 
Using our own IDL routines, we fit spectrophotometric standard star spectra to flux calibrate 
our data and remove telluric lines (I Wade fc Hornelll988l : iMatheson et al.ll2000l ). 



We give redshifts for the counterparts and nearby companions in Table |H In the cases of 
RTs 19840613, 19870422, 19920826, and 19970528, the redshift is for the galaxy that the RT 
is situated in or closest to. RT 19860115 has two possible hosts in the Keck image; Objects 
1 and 2 are located to the northeast and southwest of the RT, respectively. 



-'^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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2.4. PAIRITEL Observations 



Data were obtained witli PAIRITEL (IBIoom et al.ll2006l ) on 16 Feb. 2006. Images were 



taken simultaneously in the J, if, and Ks bands in four different pointings with a total 
integration time of about 1120 s per pointing. Additionally, to get a deep infrared (IR) 
image in one portion of the radio field, several hours of imaging in June and July 2006 were 
obtained centered at J2000 position a = 15''02"^36.0^ 6 = +78°16'39". The approximate 
3a upper limits in the shallow, wide-angle map are J = 19.2 mag, H = 18.5 mag, and 
Kg = 18.0 mag (the limits are significantly worse toward the edges of the mosaic). The 
pointings were stitched together using SWARP after finding preliminary WCS solutions in 
all bands. The result is a 1050 x 1050 pixel set of mosaic images with 1" x 1" resolution that 
overlaps most of the radio image. An astrometric solution for the mosaic images was found 
through a cross-correlation with 43 stars in common with the USNO Bl.O catalog. Typical 
rms uncertainty of the astrometric tie to the ICRS is 250 mas. 



3. Results 

We have detected a total of 8 transient radio sources in single epochs of this survey 
(Figure H]). Additionally, we have detected two transients in two-month averages of the data. 
Our source identification methods described previously argue that no more than one of these 
sources is a false detection. All of these sources are adequately fit as point sources. None 
of the transients is detected in epochs other than the detection epoch, or in longer-term 
averages with the possible exception of RT 19840613, in which the host galaxy is detected 
in the full 20-year image. 

In Figures [5] and [6] we show images of the transients overlaid on the deep radio image, 
on multi-color images from the Keck imaging, and, where Keck data are missing for two 
objects, on Kg-hand images from PAIRITEL. 

We identify several hosts or potential hosts to the RTs using the Minnesota Palomar 



Plate Survey (MAPS, http://aps.umn.edu iPennington et al.lll993l ). MAPS is a catalog of 



sources found with and characterized by plate- measuring equipment surveying the Palomar 
Sky Survey I (POSS I). MAPS coordinates are consistent with the radio reference frame 
with an accuracy of 0.3". We confirmed this by comparison of the position of two bright 
reference stars in the MAPS catalog with positions in the Tycho catalog. Referring the 
Tycho positions to epoch 1950 (appropriate for POSS), we find agreement of ~ 0.1" for GSC 
04562-00668 and ~ 1" for GSC 04562-00458. 



- 10 - 



3.1. A Transient Associated with a Spiral Galaxy: RT 19840613 



One transient is of particular interest because it is coincident with a spiral galaxy at 
a redshift z = 0.040 (a distance of ~ 170 Mpc using Hq = 72km~^Mpc~^). The galaxy, 
identified as MAPS-P023-0189928, has R = 16.5 mag, B — R = 1.8 mag, and a major axis of 
~ 20". The transient is clearly non-nuclear; it is offset to the southwest of the optical nucleus 
by 2.7" ± 0.4", corresponding to a projected distance of 2.3 ± 0.3 kpc. The optical nucleus 
is coincident with the infrared nucleus as identified in the 2MASS catalog at an accuracy of 
0.2". The probability of a single chance association with an optical source in the POSS at a 
few arcsecond precision is ^ 1%. We calculate an isotropic radio luminosity of 2 x 10^^ erg 
s~^ Hz~^, or 2 X 10^*^ erg s"^ for a bandwidth of 10 GHz. 

The transient is not detected in any other individual epoch, even with reduced detec- 
tion thresholds. The transient is also not detected in any two-month average in the year 
surrounding its detection, nor is it detected in any annual average from the entire data set. 
It is also not detected in the year prior to detection, nor in the six years following detection. 
The absence of a detection in other light curves at 5 GHz is consistent with a rapid rise and 
a decline in the flux density with a power law of a > 0.5. 

The deep image of the region surrounding the transient shows extended radio continuum 
emission associated with the spiral galaxy and a companion galaxy (Figure [7]). The location 
of the transient coincides with a region of peak emission, but it is not possible to conclude 
whether that emission is the result of the transient or a peak in the galactic emission. The 
peak of emission near the position of the transient is ^ 28 /iJy, which corresponds to a 



star- formation rate of < O.IMq yr ^ (e.g. jYun fc Carillill2002l ) 



The association of the transient with the non-nuclear region of a spiral galaxy, together 
with its luminosity, suggest that the transient may be a RSN. The transient is very unlikely 
to be a Type 11 RSN, as these sources exhibit rise and decay timescales of month s to years. 



imply ing that a Type 11 RSN would be detected multiple times in ou r survev (lEck et al. 



2OO2I ). Type la RSN have never been detected at radio wavelengths (IWeiler et al.ll2004l ). 
Type Ib/c RSN, on the other hand, evolve on timescales of days to months. The anomalous 
Type Ib/c RSN associated with SN 1998bw and GRB 980425 do ubled its flux density during 
its rise and halved its flux density following the peak in < 10 d (IKulkarni et al.lll998l ). The 
light curve declined as with the exception of a brief period of increase attributed to 
dumpiness in the circumstellar medium. The luminosity of RT 19840613 pla ces it among 



the b rightest SNe Ib/c, although a factor of a few fainter than SN 1998bw (IBerger et al. 



2OO3I ). The luminosi ty of SN 2002ap, a lso a SN Ib/c, is four orders of magnitude less than 
that of SN 1998bw JSereer et aDl2003h . The luminosity of RT 19840613 is three orders of 
magnitude less than the typical radio luminosity of a GRB afterglow. 
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3.2. A Transient Associated with a Blue Galaxy: RT 19870422 

RT 19870422 is located within ~ 1.5" of the centroid of MAPS-P023-0189613. The Keck 
image shows that this RT is clearly associated with a blue galaxy and is significantly offset 
from the nucleus. The integrated R magnitude of the galaxy is 20.2, with B — R = 2.5 mag, 
and the redshift is determined to be 2; = 0.249 (or a distance of 1050 Mpc). Galaxy-template 
matching to the Keck spectrum shows it to most likely be an Sc galaxy with strong [O II] 
and [O III] emission, indicating current star formation. 

The inferred isotropic radio luminosity is Li, 7 x 10^^ erg s^^ Hz^^, 01 L ^ 7 x 
^q39 gj^g g-i ^ bandwidth of 10 GHz. RT 19870422 was detected in a two-month inte- 
gration, implying a total energy release E ^ A x 10^^ erg. We see no evidence for detection 
in other epochs, implying a source that fades with a timescale of two months or less follow- 
ing its peak. In the deep integration, we see that the galaxy is detected with an irregular 
morphology (Figure [S])- The faint radio emission, however, is offset from the position of RT 
19870422. 

RT 19870422 shares many characteristics with RT 19840613. The RT position is clearly 
separated from the peak of the optical emission, ruling out the possibility of an AGN event. 
Its longer duration and higher luminosity, however, differentiate it. The luminosity is nearly 
an order of magnitude greater than that of S N 1998bw, and two orders of magnitude more 



luminous than the bright Type 11 SN 1979C (lEck et al.ll2002l ). Thus, the luminosity falls in 



between the maximum observed luminosity of a RSN and the typical luminosity of a GRB 
afterglow. The two-month duration of RT 19870422 is more consistent with known RSNe of 
both Type Ib/c or Type 11, as well as GRB afterglows. 



3.3. Transients Possibly Associated with Galaxies or Galaxy Groups: RT 

19920826 & 19970528 

RT 19920826 is located ~ 5" to the southeast of a red galaxy identified in the Keck 
image, which is not identified in the MAPS catalog. The host is a red object with no clear 
emission or absorption lines. 

RT 19970528 is located ~ 5" to the northwest of MAPS-P023-0189499, which has R = 
19.6 mag and B — R = 2.4 mag. Two other galaxies are identified within a radius of 30". 
The nearest one is a red, early-type galaxy at z = 0.245. Galaxy-template matching shows 
it is most likely an elliptical galaxy, but it is also well fit by an SO spectrum. 

We computed the probability of false association of these galaxies with transients through 
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Monte Carlo simulations. We randomly distributed 100,000 radio sources over the field and 

then measured the distance to the nearest galaxy as identified in the MAPS catalog. The 
probability of a false association at a distance of less than 5" is only 2%. The probability of 
a false association at a distance of 15" is ~ 20%. The probability for any one of the 10 RTs 
being randomly associated with a galaxy at a distance of 5" is ~ 20%, and at a distance 
of 15" the probability is nearly unity. The probability of multiple sources located within 5" 
is ^ 5%. Thus, the association of either RT 19970528 or RT 19920826 or both RTs with a 
galaxy is probable but not certain. 

3.4. Sources without Radio or Optical Counterparts: RT 19840502, 19860115, 

19860122, 19970205, 19990504, & 20010331 

There is no evidence for faint radio, optical, or infrared counterparts for six of our 
sources: RT 19840502, 19860115, 19860122, 19970205, 19990504, 20010331. There are no 
sources apparent in the image obtained by averaging the succeeding two months of data, as 
well as in the year-long and decade-long averages. These results require a rapidly decreasing 
light curve, proportional to or steeper. 

RT 19860115 has two possible hosts. The first possible host. Object 1, is located ~ 15" 
to the northeast of the RT and is a starburst galaxy a,t z — 0.1297 with strong emission 
fines and a blue continuum. Object 1 is identified as MAPS-P023-0190130, a galaxy with 
R = 19.2 mag and B — R = 1.3 mag, and is also detected as a faint, extended radio source. 
The other possible host. Object 2, is located ~ 20" to the southwest of the RT and is an 
early-type galaxy ai z = 0.242. Galaxy-template matching shows the spectrum is best fit by 
either an elliptical or SO galaxy spectrum. Several other objects, likely to be galaxies, are 
in the field. The physical separation of RT 19860115 from the two objects is ~ 50-100 kpc, 
which is still within a radius at which the galaxy may contain stars. Nevertheless, given the 
Monte Carlo simulations discussed above, we conclude that this transient is likely associated 
only by chance with these galaxies. 

For the remainder of sources, there are no obvious associations with galaxies or stars. 
The absence of optical or infrared counterparts only eliminates host candidates; the non- 
simultaneity of radio and optical observations implies that we place no limits on optical 
transients. 
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4. Radio Transient Rates 

We can compute rates of transients in several different ways. We are limited in the 
accuracy of estimating transient rates by our lack of knowledge of the characteristic time 
scale, tchar, of the transient sources. This characteristic time falls between the duration of 
the observations (~ 20 min) and the typical separation between epochs (7 days). 

The two-epoch survey sensitivity is independent of tchar- The effective area for a two- 
epoch survey for a given flux-density limit is just the sum of the areas in each epoch in which 
a source above the flux-density limit is detectable. For a set of images with uniform 
sensitivity over a solid angle Q in which Nt transients are detected, the two-epoch source 
density is Nt{Nf. — l)~^n~^. For a sample with varying image properties, we can compute 
the two-epoch source density given the image noise statistics, the synthesized beam size, 
and the primary beam shape. The ratio of the primary beam solid angle to the synthesized 
beam solid angle gives the number of independent pixels in the image. The source detection 
threshold per image is then set based on the number of independent pixels and the rms noise. 
We combine the 5 GHz and 8.4 GHz survey limits and transient detections together. The 
accumulated results as a function of fiux-density threshold are listed in Table [51 



We use a non-parametric Kaplan-Meier method (IFeigelson fc Nelsonlll985l ). for estimat- 
ing the cumulative source-count distribution and then compute transient rates (Figure [H]). 
Seven transients are detected above a flux density of 370 yuJy, giving a two-epoch source 
density of 1.5 ± 0.4 deg~^. The transient rates are reasonably fit by a distribution with S~'^ 
and 7 = 1.5, which corresponds to a non-evolving, volume-limited sample. However, the 
constraints on 7 are not very strong. A larger 7 appears to be required to ac hieve agreement 



between the rate determined from this survey and the NVSS-FIRST survey (iGal-Yam et al. 



20061 ). Different event rates at 1.4 and 5 GHz due to synchrotron self-absorption or other 
spectral index effects may also influence this comparison. Additionally, the FIRST-NVSS 
survey may be incomplete due to confusion from the mismatched resolution of these surveys. 

We detected no transients in the 17 annual average images. The typical flux detection 
threshold at the half-power radius is 90 yuJy. The total two-epoch area surveyed at this 
sensitivity is 0.3 deg^, implying a 2cr upper limit to the two-epoch density for year-long 
transients at 90 /iJy of ~ 6 deg~^. In 96 images composed of two- month averages, we found 
two transients. At a flux-density threshold of 200 /iJy, the effective area surveyed is 1.9 deg^, 
giving a 2a upper limit to the two-month transient rate of ~ 2 deg~^. 

We can also determine a rate as a function of area and time. The precision of this rate 
is limited by our lack of measurement of the characteristic time of the transient events. The 
rate for sources brighter than 370 /iJy is then 0.07 < i? < 40 deg~^yr~^. 
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If we assume that RTs 19840613 and 19870422 belong to a distinct transient class based 
on their association with blue galaxies, we can estimate a rate for such transients as a function 
of volume and time. We could detect this transient to a distance of ~ 240 Mpc, which gives 
a total volume of ~ 300 Mpc^. Assuming tchar = 7 d, we find a rate of 3 x 10^ Gpc~^yr~^. 
Similarly for RT 19870422, we find a maximum volume of ~ 5000 Mpc^ and use tchar = 60 d 
to find a rate of 4 x 10^ Gpc-^yr"^- Jointly for RT 19840613 and RT 19870422, the rate is 
~ 2 ± 2 X 10^ Gpc^^yr^^. These rates are consistent wi thin an order of rn agnitude of the 
SN Ib/c rate determined optically, 4.8 x 10^ Gpc~^yr~^ (IBerger et al.ll2003l ). 



There is some evidence that our VLA field contains an overdensity of galaxies at z ~ 0.25 
(Tabl e HI). Additionally , the cluster Abell 2047 is located within 9' of the center of the VLA 
field ( lAbell et al.lll989l ). A2047 is a low-density cluster (richness class 0) with a radius ~ 5' 
and a distance class of 6. The distance class is an estimate of distance based on visual 
inspection of photog raphic plates. The di stance class of 6 is consistent with a redshift of 
0.25 for the cluster (IStruble fc RoodI 119871 ). If all the RTs are extragalactic, then the RT 
rates may be overestimated relative to the rate for the field by no more than a factor of a 
few. 



5. Discussion 



Two of the transients are reasonably identified as RSNe, possibly Type Ib/c or II and 
similar to known RSNe associated with GRBs. We consider now whether the eight transients 
not clearly identified with nearby galaxies are examples of known or expected classes of radio 
transients. These transients are unlikely to be RSNe unless they are associated with faint 
and/or low surface brightness galaxies. As we demonstrate below, the observed transient 
parameters do not easily fit into any expected class. 

The known and anticipated sources of radio transients are very broad. Since we have 
no evidence for variability on a timescale of less than 20 min, we exclude from further 
discussion transient sources that have characteristic timescales of les s than 1 s. These include 
pulsars and related phenomena such as RRATS and giant pulses (IMcLaughlin et al.l 12006 



Hankins et al 



(iBastian et al. 



2003 



2000 



Lazio et al 



200 



IT 



Cordes fc McLaughlinI |2003| ). as well as flares from extrasolar planets 
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5.1. X-ray Sources 



Information on the X-ray properties of the transients and transient hosts is very hmited. 
The ROSAT catal og of bright sources does not identify any sources in the field of our survey 
(IVoges et al.lll999l ). There are also no sources in the RXTE All Sky Monitor catalog in the 
field, in dicating an u pper limit to the X-ray flux of events between 1997 and 2001 of ~ 10 
mCrab (jSwankI Il999l ). No observations of the field have been made with the Chandra or 
XMM Newton X-ray Observatories. 



5.2. Gamma-Ray Burst Afterglows 



Higher-energy counterparts to these transients have not been clearly detected. There 
are no known optically discovered SNe in the field of view of our survey, nor have there 
been any optical variability campaigns in this field. There are no GRBs which are uniquely 
localized to this field, although there are bursts with large positional uncertainties that could 
have originated in the field of view of our survey. We identify 7 BATSE events that occurred 
within la of our field center between 1991 and 1999. Uncertainties in the BATSE positions, 
however, are ~ 10°, indicating a probability of < 0.1% for even one of these GRBs to have 
occurred within our field of view. Of the GRBs, one event (GRB 990508) occurred 4 d 
after RT 19990504. Given the low probability of concurrence and the lack (in other GRBs) 
of radio emission before gamma-ray emission, we conclude that this is almost certainly a 
coincidence. There is no systematic relationship between the remainder of the events and 
the radio transients. 



5.3. Orphan Gamma-Ray Burst Afterglows 



In the standard paradigm of GRBs and GRB afterglows, the gamma-ray emission is 
highly anisotropic while the afterglow emission is less anisotropic or even isotropic. Thus, 
the observed GRBs and their afterglows from Earth are a subset of the total number of 
GRBs. Afterglows detected without gamma-ray emission are known as orphan gamma-ray 
burst aft erglows. The rati o of OGRBAs to GRBs is the inverse of the beaming fraction. 



/{, ^ (e.g.. lDalal et al.ll2002l ). The phenomenon of a GRB afterglow without a corresponding 



gamma-ray burst has not been conclusively detected. 

Detailed models of OGRBAs predict ~1 per square degree at se nsitivity thresholds of 



0.1 mJy, which is consisten t with our number of transient detections (ILevinson et al.l 12002 



Totani &: Panaitescul |2002| ). Whereas the canonical radio afterglow for the forward shock 
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rises on a timescale of days and decays slowly, the majority of the radio transients found in 
the present survey occur on timescales less than 1 week. One possible explanation in the 
context of GRB orphans is that the emission is due to the reverse shock, which is known, 
in the few cases where follow-up observati ons occurred soon after the GRB , to produce a 
bright and r apid (< 3 d) rad io signature (iKulkarni et al.l Il999l : iFraill l2003l . e.g.,). Strong 
scintillation (IFrail et al.lll997l . e.g.,) of a standard, faint afterglow could also be responsible 
for producing short-term amplification that leads to an apparent fast transient, but then we 
might see marginal evidence for the forward shock in the epochs following the transient — 
yet this is not seen in individual sources. 

The single long-timescale transient (RT 20010331) could be consistent with having arisen 
from the long-lived for ward shock (as it tran sitions from mildly relativistic to nonrelativistic). 
Levinson et al. (2002) iLevinson et al.l (120021 ) calculate the expected number of forward-shock 
dominated OGRBAs for radio surveys. Taking the canonical interstellar medium density of 
no = 1 cm~^, it is reasonable to expect that we only find 1 viable forward-shock candidate 



with the nominal beaming of 



500. Nevertheless, since for the parameters of our 



survey we probe only to 2 ~ 0.6, the non-detection of an apparent host galaxy in deep Keck 
imaging of RT 20010331 would seem to suggest that our basic model assumptions for the 
radio afterglow are suspect. 

Thus, the transients that we have detected can be considered consistent with the ex- 
pectations of OGRBAs. There is, however, a large diversity of high-energy phenomena 
with poorly explored electromagnetic signatures beyond the classical G RB afterglow pic- 



ture, including X-ray flashes , X-ra y rich GRBs (ISakamoto et a. 



2005). GRBs associated 



with SNe Ib/c (iBerger et al.ll2003D . and short-duration GRBs fiGehrels et al.l 120051 1. The 
richness of this phenomenology suggests that the radio transients we have discovered could 
reflect the continuum of stellar collapse events. 



5.4. Active Galactic Nuclei 



AGNs are highly variable and common. With the appropriate selection techniqu e, AGNs 
can be found with a number density of > 100 per square degree (jStern et al.l 120051 ). If our 
transients are AGNs, then the absence of deep, static detections indicates they are associated 
with radio-quiet objects, which constitute the majority of AGNs. Short-timescale variability 
of AGNs is often associated with inter stellar scintillation and extreme scattering events, 
and typically has amplitudes of ~ 10% (iLovell et al.ll2003l ). Variations of a factor of a few 
on timescales of hours have been seen in the most extreme cases such as PKS 0405-385 and 
J1819+3845, which is not sufficient to account for the observed transients given the very large 
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amplitude of variation observed (IRickett et al.ll2002l : iDennett-Thorpe fc de Bruynll2000l ) 



Intrinsic variations am ong luminous AGN s are observed to have a timescale of months to 
years at radio frequencies (iHughes et al.lll992l ). Variations are rarely larger than a factor of 
a few in radio-loud objects. A survey of variable background sources in the direction of M31 
found a small number of o bjects with variability of a factor of a few to 10 on timescales of 
years (IGelfand et al.ll2005l ). An interesting example is the nucleus of the nearby spiral galaxy 
III Zw 2, which has exhibited factor of 20 fluctuations in its radio luminosity on timescales 



2005h The T eV blazar Mrk 421 shows 50% variations 



of months to years (IBrunthaler et al. 
on a timescale of ~ 10 d (IBlazejowski et al.ll2005l ). Sudden increases in the accretion rate, 
or shocks in the jet, could lead to a dramatic flare at radio wavelengths; however, these 
flares are expected to have timescales on the order of months to a year, inconsistent with 
our observations. 

Recently discovered evidence for tidal disruption events in AGNs observed in the u l- 
traviolet indicate the possibility for large-amplitude AGN transients (IGezari et al.l 120061 ). 
However, the timescale for the observed UV transients are months to a year rather than 
days. Without a mechanism for a short timescale in the radio, tidal disruption events are 
unlikely to account for the observed radio transients. 



5.5. Stellar Sources 

We cannot exclude a Galactic origin for these transients. A wide range of stars and 
stellar systems are known to flare on timescales of minutes to hours to days. However, no 
known stellar types closely fit the properties of these transients. Our field (6 = 36°, / = 115°) 
is out of the plane of the Galaxy, away from the Galactic center and bulge, and is not near 
any known star cluster or star-forming region. It is likely, therefore, that any stellar origin 
for these transients would be even more numerous in one of these more favored directions. 



5.5.1. Known Luminous Transients 



The luminosity of the observed transients is Ly 10^^((i/l kpc)^ erg s ^ Hz ^. T his is 



near the upper limit of the known stellar radio luminosity distribution (jGudelll2002l ). RS 
CVn binaries, FK Com class st ar s , Algol- class st ars , and T Tauri stars have been observed 



to radiate at this luminosity (e.g.. lBower et al.ll2003l ). The total number of known RS CVn 



FK Com stars, and Algol stars discovered at X-ray wavelengths is < 1000, indicating that we 
are unlikely to detect any such objects by chance, even if the distance cutoff in our survey is 
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10 kpc (jPempsey et al.lll993l ). T Tauri stars are more common but are unlikely to be found 
outside regions of active star formation. 



5.5.2. Late-Type Stars 

M-type dwarfs at distances as lar ge as 1 kpc have been identified as an i mportant 
contribution to optical transient rates (IBecker et al.ll2004l : iKulkarni fc Raull2006l ). Due to 



the typical low radio luminosity {Ly ^ 10 - 10 erg s~ Hz~ ) of known dMe stars and 
brown dwarfs, they are likely to be detected at distances of 30 pc or less. In an extreme case, 
the dMe star EV Lac was observed to undergo a flare of a factor of ~ 150 to a maximum 
l uminosity L,, k, 1 0^^ erg s~^ Hz~^, with a rise time of minutes and a decay time of hours 



( lOsten et al.l 120061 ). Such a flare would be observed in our survey to a distance of 100 pc. 
The local space density of stars with M < 0.3 Mq is 0.035 pc~^, implying only a s mall chance 



of a low-mass star or brown dwarf within the volume to which we are sensitive (IReid et al. 



19991 ). Any stellar origin for these transients must come from a faint source, since we have 
no detection of stellar counterparts in Keck, 2MASS, or PAIRITEL images. The 2MASS 
magnitude limit of Kg = 15 implies that any L dwarf {Mk ~ 10 mag) or earlier-type star 
must be at a distance > 100 pc to remain undetected. PAIRITEL limits require such stars 
to be at a distance > 500 pc. The limiting Keck magnitude requires late-type M stars to be 
at a distance > 1000 pc. 

The probability distribution of radio activity in low-mass stars is not well quantified. 
Few stars have been observed at radio wavelengths for more than a few hours or a few days. 
Despite the lack of long-term monitoring of dMe stars, one can estimate the probability 
distribution of flares by comparison with the properties of the Sun. The prob ability of a 
solar flare at centimeter wavelengths of a given flux density S scalesasjS"^''^ (INita et al. 



2004 ). /^From observations of a large sample of late-type stars (lBergerll2006l ). we estimate 
a probability of ~ 0.01 for a late-type star to exhibit a luminosity L^, 10^^ erg s~^ Hz~^. 
If the probability of stellar flares scales with the same relation as solar flares, a flare of 
Ly ^ 10^^ erg s~^ Hz"^ from a late-type star at a distance of 1 kpc would be visible in our 
radio survey but missing in the Keck image, and it would occur with a probability of < 10^^. 
Assuming a constant density of stars set to the local value, the total number of M dwarfs 
out to 1 kpc within our field of view is < 5 x 10'^. Given our ~ 1000 observations, the 
expected number of flares of this magnitude detected is ^1. This is an upper limit since 
the statistics for the observed stars are biased by selection for activity at other wavelengths. 
Within an order of magnitude, though, this matches the number of transients that we detect, 
suggesting that low-mass stars at 1 kpc or greater distance may contribute to, but are unlikely 
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to dominate, the observed radio-transient rate. 



Very large radio flares might be accompanied by luminous X-ray flares. If the transients 
follow the radio/X-ray correlation (jGudeJ |2002) , then a radio flare of 10^*^ erg s~^ Hz~^ at 
a distance of 1 kpc would have an X-ray luminosity of 10~^^^ erg s~^, which exceeds the 
RXTE ASM detection threshold of ~ 10~^° erg s~^. Howeve r, many sources do not closely 
follow this relationship, including many low-mass stars (e.g.. iBerger et al.ll200ll ). 



5.5.3. Soft Gamma-Ray Repeaters 
The transients share some properties with soft gamma-ray repeaters (SGRs) . SGRs are 



magnetars that exhibit dra matic flares at gamma-ray and radio wavelengths (e.g.. lFrail et al. 



19991 : iGaensler et al.ll2005l ). SGR 1806-20 produced a flare that reached a luminosity of 
5 X 10^^ erg s~^ Hz~^, had a characteristic timescale of days, and declined as t"^'''. Such a 
flare would be detectable in our dataset at a distance of ~ 200 kpc, would appear in only 
a single epoch, would fade away rapidly so that it is not seen in long-term averages, and 
would have no lasting optical counterpart. If all of our transients are Milky Way SGRs to 
a distance of ^ 200 kpc, the implied rate is ^ 0.01 kpc~^ y~^ at high Galactic latitude. 
Without making any correction for Galactic population distribution, we estimate a total 
rate of ^ 3 yr~^ in the disk of the Galaxy, which is at least an order of magnitude too high 
relative to the number of discovered SGRs. 



5.5.4. X-ray Binaries 



X-r a.y binaries (XRBs ) can be seen throughout the Galaxy (e.g.. lBower et al.ll2005l ) but 
are rare (IMuno et al.ll2005l ). For the extremely overdense Galactic center, the XRB detection 
rate in the radio is ^ 0.01 kpc^'^ yr^^. The comparable observed transient rate would imply 
a similar overdensity at high Galactic latitude, which is implausible. 



5. 5. 5. Pulsars 



Most pulsar variability occurs on time scales much shorter than our observations and is 
therefore unlikely to account for our observed transients. There have been some examples, 
however, of longer timescale variability. PSR B1931-I-24, for insta nce, has periods of a ctivity 
and inactivity with characteristic times of days to tens of days (IKramer et al.ll2006l ). This 
variability is quasi-periodic, however, indicating that we would detect it and similar pulsars 
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in multiple epochs throughout our survey. 



5.6. Propagation Effects 

We have already referred to the phenomenon of intra-day variability (IDV), which in 
the most extreme cases leads to variations of a factors of a few on a timescale of hours. IDV 
of this scale has been demonstrated to orig inate from interstellar scintillation of compact 



structure in the source. iGalama et al.l (119971 ) identified a factor of 43 fluctuation in the flux 
density of PSR B0655+64 at 325 MHz, which they argue originates from strong focusing in 
the interstellar medium — that is, a caustic. It is difficult to estimate the frequency with 
which such a phenomenon can occur, but observational evidence is scarce. 

Microlensing could also produce a transient amplification of a background source. A 
simple calculation estimates that to explain all 10 RTs as the result of microlensing requires 
the product nr ^ 10~^, where n is the number of sources that can be lensed and r is the 



lensing opacity (jPaczynskil Il996l ). Since r has been measured to be ~ 10~^ by MACHO 
and OGLE experiments, then n ^ 10^, which is absurdly large for any known population of 
Galactic or extragalactic radio sources. 

Finally, refiected solar fiares off objects in the Solar System could produce detectable 
radio flares. A 1 MJy solar flare reflecting off a 1000 km object at a distance of 1 AU from 
the Earth would produce a 0.1 mJy flare at Earth. Our observed fleld is far out of the 
ecliptic plane, however, and there are few known Solar-System objects of this size outside of 
the ecliptic. Nearer objects could be smaller but would move significantly during the course 
of the observation. 



6. Conclusions 

We have conducted a 944-epoch, 20-yr-long survey for radio transients with archival 5 
and 8.4 GHz Very Large Array data. Ten radio transients are apparent in this data set. 
Eight transients appear in only a single epoch and disappear completely thereafter. Two 
transients appear in two-month integrations and also never reappear. We estimate the rate of 
radio transients below 1 mJy with a large uncertainty due to the wide range of characteristic 
times possible for these sources. 

Two of the transients may be peculiar Type Ib/c and/or Type II radio SNe in nearby 
galaxies. All of these transients share some characteristics with those expected of orphan 
gamma-ray burst afterglows, but they cannot be conclusively identified as such. The diversity 
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of high-energy stellar-death phenomena makes it difficult to place specific limits on the 
nature of GRBs based on these results. The absence of persistent counterparts to these 
sources indicate that they are unlikely to be AGNs. Galactic sources appear unlikely to 
be an explanation, but they certainly cannot be ruled out. In particular, late-type stars 
could contribute to the number of radio transients if their luminosity distribution can be 
extrapolated by several orders of magnitude beyond the most luminous observed M-dwarf 
flare. Deeper multi-wavelength imaging is critical for identifying the hosts of these events. 



In the near future, the Allen Telescope Array (iDeBoer et al.l 12004 ) will be efficient at 
discovering transients of this sort. The wide field of view of this telescope makes it a powerful 
instrument for large solid angle radio surveys. Given the apparently short timescale of these 
phenomena, rapid response to transient discovery at radio and other wavelengths will be 
critical for determining the nature of these transients. 
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Table 1. Steady Sources Detected in at Least One Epoch 



RA 




Dec. 




Ndet 


< 5> 


c 

^deep 


(J2000) 




(J2000) 




(^Jy) 


(MJy) 








5 GHz 








15 00 11.48 ± 0.10 


78 


12 12.26 


± 00.13 


26 


2380 ± 73 


502 ± 17 


15 00 44.87 ± 0.10 


78 


18 39.04 


± 00.32 


4 


874 ± 36 


180 ± 7 


15 01 11.70 ± 0.10 


78 


15 20.23 


± 00.10 


141 


504 ± 6 


307 ± 4 


15 01 16.19 ± 0.10 


78 


12 45.88 


± 00.10 


363 


1146 ± 14 


741 ± 5 


15 01 22.67 ± 0.10 


78 


18 05.66 


± 00.10 


452 


880 ± 9 


634 ± 3 


15 02 51.40 ± 0.10 


78 


18 12.46 


± 00.10 


23 


351 ± 19 


155 ± 3 


15 03 24.83 it 0.14 


78 


17 37.62 


± 00.46 


1 


318 ± 54 


166 ± 4 


15 03 28.07 ± 0.30 


78 


09 21.50 


± 01.04 


1 


2045 ± 301 


2650 ± 8 








8.4 GHz 








15 01 11.13 ± 0.12 


78 


15 22.45 


± 00.32 


1 


1081 ± 154 


291 ± 9 


15 01 22.58 ± 0.10 


78 


18 05.80 


± 00.10 


159 


1099 ± 36 


794 ± 8 


15 02 51.38 ± 0.10 


78 


18 13.14 


± 00.18 


3 


518 ± 31 


225 ± 4 


15 03 23.41 ± 0.41 


78 


17 35.74 


± 00.64 


1 


840 ± 157 


212 ± 8 



Note. — Columns arc (1) right ascension, (2) declination, (3) the number of epochs in which the source is detected, (4) mean 
flux density of detections, and (5) flux density in the deep image. 



Table 2. Transient Sources from a Single Epoch 



Epoch 


RA 


Dec. 


S 




tnext 


Snext 


-i( 


agio Host 




(J2000) 


(J2000) 


(MJy) 


{/"Jy) 


(days) 


(MJy) 













5 GHz 










19840502 


15 02 24.61 ± 0.21 


78 16 10.08 ± 00.46 


448 ± 74 


< 8 


7 


-10 ±68 


4.5 


19840613 


15 01 38.07 ± 0.21 


78 18 40.75 ± 00.39 


566 ± 81 


< 28 


7 


86 ±89 


7.5 


19860115 


15 02 26.40 ± 0.44 


78 17 32.39 ± 01.59 


370 ± 67 


< 8 


7 


199 ± 121 


3.6 


19860122 


15 00 50.15 ± 0.34 


78 15 39.37 ± 01.38 


1586 ± 248 


< 15 


7 


-59 ± 164 


6.4 


19920826 


15 02 59.89 ± 0.35 


78 16 10.82 ± 02.44 


642 ± 101 


< 9 


56 


37 ±83 


6.3 


19970528 


15 00 23.55 ± 0.10 


78 13 01.37 ± 00.17 


1731 ± 232 


< 36 


7 


90 ± 206 


7.9 


19990504 


14 59 46.42 ± 0.56 


78 20 29.03 ± 00.74 


7042 ± 963 


< 117 


21 


-313 ± 1020 


9.1 








8.4 GHz 










19970205 


15 01 29.35 ± 0.10 


78 19 49.20 ± 00.10 


2234 ± 288 


< 646 


5 


857 ± 323 


8.0 



X 
G 
X 
X 

? 
? 
X 



Note. — Columns are (1) RT epoch, (2) right ascension, (3) declination, (4) mean flux density of detection, (5) flux density 
in the deep image, (6) separation in days between next epoch and epoch of transient, (7) flux density in the next epoch, (8) the 
negative logarithm of the probability of false detection in the single epoch, and (9) optical/IR host identification, where "G" 
indicates galaxy, "?" indicates possible galaxy association, and "X" indicates no host. 
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Table 3. Transient Sources with Tchar = 2 months 



Epoch 


RA 


Dec. 


S 


^deep 


tnext 


Snext 


- logio PFD 


Host 




(J2000) 


(J2000) 




{/"Jy) 


(days) 














5 GHz 












19870422 


15 00 50.01 ± 0.35 


78 09 45.49 ± 0.97 


505 ± 83 


< 36 


96 


-154 ± 65 


3.0 


G 








B.4 GHz 












20010331 


15 03 46.18 ± 0.10 


78 15 41.68 ± 00.11 


697 ± 94 


< 37 


59 


85 ± 85 


3.4 


X 



Note. — Columns are (1) RT epoch, (2) right ascension, (3) dcchiiatioii, (4) mean flux density of detection, (5) flux density 
in the deep image, (6) separation in days between next epoch and epoch of transient, (7) flux density in the next epoch, (8) the 
negative logarithm of the probability of false detection in the single epoch, and (9) optical/IR host identification, where "G" 
indicates galaxy, "?" indicates possible galaxy association, and "X" indicates no host. 



Table 4. Magnitudes, Colors, and Redshifts of Transient Host Galaxies 



RT Epoch 


MAPS Counterpart 


Host Offset 


R 


B-R 


z 








(mag) 


(mag) 




19840613 


P023-0189928 


2.7" NW 


16.5 


1.8 


0.040 


19860115-1 


P023-0190130 


15" NE 


19.2 


1.3 


0.130 


19860115-2 




20" SW 






0.242 


19870422 


P023-0189613 


1.5" SE 


20.2 


2.5 


0.249 


19920826 




5" NW 






unknown 


19970528 


P023-0189499 


5" SE 


19.6 


2.4 


0.245 



Note. — Host offset is position relative to the transient. Magnitudes and colors are from 
MAPS catalog. 
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Table 5. Effective Area Searctied for a Comparable Two-Epocli Survey 





5 GHz 




8.4 GHz 




Total 




Flux Density 


Number of Epochs 


Area 


Number of Epochs 


Area 


Number of Epochs 


Area 






{dcg2) 




(dcg2) 




(dcg2) 


70 


9 


0.16 





0.0 


9 


0.16 


100 


10 


0.24 





0.0 


10 


0.24 


140 


12 


0.34 





0.0 


12 


0.34 


200 


35 


0.50 





0.0 


35 


0.50 


280 


201 


1.42 


442 


1.45 


435 


2.19 


400 


464 


4.57 


560 


3.08 


761 


6.20 


560 


561 


9.14 


584 


4.79 


873 


11.68 


800 


596 


14.56 


591 


6.64 


910 


18.09 


1120 


606 


19.85 


591 


8.40 


924 


24.31 


1600 


607 


25.51 


598 


10.28 


925 


30.97 



Table 6. Effective Area Searched for a Comparable Two-Epoch Survey for Tchar = 2 

months 





5 GHz 




8.4 GHz 




Total 




Flux Density 


Number of Epochs 


Area 


Number of Epochs 


Area 


Number of Epochs 


Area 






(deg2) 




(deg2) 




(deg2) 


70 


24 


0.10 


9 


0.02 


26 


0.1 


100 


41 


0.49 


38 


0.22 


54 


0.56 


140 


55 


0.91 


61 


0.67 


73 


1.11 


200 


75 


1.51 


74 


1.31 


99 


1.91 


280 


90 


2.24 


78 


1.97 


115 


2.84 


400 


95 


3.11 


78 


2.69 


120 


3.94 


560 


95 


3.93 


78 


3.37 


120 


4.97 


800 


95 


4.81 


78 


4.09 


120 


6.07 


1120 


95 


5.64 


78 


4.77 


120 


7.11 


1600 


95 


6.52 


78 


5.48 


120 


8.22 
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Fig. 1.— Image from epochs 19840613 (left) and 19840620 (right) at 5 GHz. The two 
westernmost sources are steady sources J1501 12+781520 and J150123+781806. The north- 
ernmost source is a transient, J150138+781841, detectable only in epoch 19840613. Contours 
are 5, 6, and 7 times the image rms of 57 and 63 /iJy. 
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Fig. 2. — Deep image obtained at 5 GHz. We show the positions of single-epoch transient 
sources at 5 GHz (squares), single-epoch transient source at 8.4 GHz (triangle), two- month 
epoch transient sources (diamonds), and steady sources detected at least once in an individual 
epoch (circles). The two large circles indicate the half-power radius and twice the half-power 
radius of the 5-GHz field. The grey scale stretches from —2.5 to 30 //Jy. The image has not 
been corrected for primary beam attenuation. 



-33- 





Fig. 3.— Light curve of J150123+781806 at 5 GHz (top) and 8.4 GHz (bottom) from the 
entire data set. 
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Fig. 4. — Light curves of the transient sources for the year surrounding detection. A circle 
around the point denotes the epoch of detection. Other points for all sources are all non- 
detections; values given are best-fit Gaussians at the position of the transient. 
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Fig. 5. — Contour plots of transients overlaid on color images from the deep radio image. 
Each image is 1' in scale. Contours are 3, 4, 5, 6, and 7 times the image rms. The synthesized 
beam for each image is displayed in blue in the lower left-hand corner. The color scale of 
the deep radio image has a range of to 100 /iJy. 
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Fig. 6. — Optical and infrared images at the positions of transient sources. The white circle 
indicates the uncertainty in the position of the RT. Images are from Keck except for RTs 
19990504 and 20010331, which are from PAIRITEL. The color composite for the Keck images 
was made by using the g-hand image for the blue channel, R band for the red channel, and 
the geometric mean of the two images for the green channel. Keck images are 38" on a side. 
The PAIRITEL image is in the K band and is 73" on a side. The orientation and location 
of the spectroscopic slit is indicated with blue bars. North is up and east is to the left. The 
bright streak in the image of RT 19860122 is a diffraction spike from a nearby bright star. 



78 19 05 F 



00 - 



o 



18 55 — 



O 

lU 

Q 



O 
I- 
< 



o 



o 

CM 
-9 




30 - 



25 - 



20 - 



15 01 46 



44 



42 



40 38 36 

RIGHT ASCENSION (J2000) 



34 



32 



Fig. 7. — Deep image of the region centered on the transient 19840613. The small cross 
marks the position and 3a uncertainty in its position. The position is clearly offset from the 
centroid of the radio emission and from the centroid of the optical emission, marked with 
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the large cross. The extent of the large cross delineates the major and minor axes of the 
galaxy as measured in 2MASS. Contours are 3, 4, 5, 6, 7, 8, 9, and 10 times the image rms 
of 2.6 //Jy. 
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Fig. 8. — Deep image of the region centered on the transient 19870422. We mark the position 
and 3(7 uncertainty of the RT (plus sign) and the position of the optical centroid (star). The 
RT position is offset from any peak of the radio emission but is consistent with the peak of 
the optical emission. Contours are —2, 2, 3, and 4 times the image rms of 14 /iJy at the 
location of the transient. 
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Fig. 9. — The cumulative two-epoch source density for radio transients as a function of flux 
density. The solid black line shows the rate while the dotted lines show the 2a upper and 
lower bounds. The red dashed lines show S^^'^ and S'"^'"^ curves. Both are reasonable fits 
to these data, while the latter shows greater consistency with lower-frequency results. The 
arrows show 2cr upper limits for transients from this survey with a 1-year timescale, two- 
month ti mescale, and for tran sients from the comparis o n of t he 1.4 GHz NVSS and FIRST 
surveys flGal-Yam et al.l 120061 ). from the ICarilli et al.l (120031) survey (labeled C2003), and 
from the iFrail et al.l (120031 ) survey (labeled F2003). See the main text for details of these 
surveys. 
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Table 7. Effective Area Searcfied for a Comparable Two-Epoch Survey for Tchar = 1 year 





5 GHz 




8.4 GHz 




Total 




Flux Density 


Number of Epochs 


Area 


Number of Epochs 


Area 


Number of Epochs 


Area 






(deg2) 




(deg2) 




(deg2) 


50 


5 


0.02 


10 


0.02 


15 


0.04 


70 


13 


0.09 


14 


0.05 


15 


0.10 


100 


17 


0.24 


15 


0.08 


20 


0.26 


140 


17 


0.38 


15 


0.14 


20 


0.42 


200 


17 


0.53 


16 


0.18 


21 


0.57 


280 


17 


0.67 


16 


0.23 


21 


0.73 


400 


17 


0.82 


16 


0.27 


21 


0.89 


560 


17 


0.96 


16 


0.32 


21 


1.14 


800 


17 


1.09 


16 


0.35 


21 


1.19 


1120 


17 


1.15 


16 


0.36 


21 


1.24 


1600 


17 


1.16 


16 


0.36 


21 


1.25 



